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Abstract 



New technique of integration of certain types of partial differen- 
tial equations is developed. For this purpose non-commutative inte- 
gration over Cayley-Dickson algebras is used. Applications to non- 
linear vector partial differential equations of Korteweg-de-Vries and 
Oh ■ Kadomtzev-Petviashvili types and describing non-isothermal flows of 

■ incompressible Newtonian liquids are given. 
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; 1 Introduction. 

oo 

This article is devoted to a method of non-commutative integration of sys- 
tems of partial differential equations or a partial differential equation writ- 
en ; 

ten in a vector form over Cayley-Dickson algebras. This technique is applied 
to non-linear generalized Korteweg-de-Vries and Kadomtzev-Petviashvili and 
also describing non-isothermal flows of incompressible Newtonian liquids par- 
^ \ tial differential equations written in vector forms with large number of vari- 

ables. Dirac had used complexified quaternions to solve Klein-Gordon's hy- 
perbolic partial differential equation, which is used in quantum mechanics. 
This approach is generalized in this paper for non-linear partial differential 
equations over the Cayley-Dickson algebras. 

The method based on R-linear integral equations was previously used for 
non-linear partial differential equations with a real time variable t and real 
space variables x, y [TJ [T7]. Using the non-commutative line integration over 
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the Cayley-Dickson algebras it is spread in this paper for variables x, y in the 
Euclidean space R", n G N. 

Korteweg-de-Vries equations describe a motion of shallow waters, for ex- 
ample, in channels. Kadomtzev-Petviashvili and modified Korteweg-de-Vries 
equations are used in modeling of physical processes such as non- stationary 
spread of waves in a material with dispersion, magneto-hydrodynamical waves 
in a non-collision plasma, heat conductivity of a lattice of anharmonic os- 
cillators. Moreover, partial differential equations describing stationary non- 
isothermal flows of incompressible Newtonian liquids are very important as 
well. This method is also applicable to a resonance interactions in a non- 
linear material, for example, in the non-linear optics. 

In this article previous results of the author on functions of Cayley- 
Dickson variable and non- commutative line integrals over Cayley-Dickson 
algebras are used [HI EJ QUI EE] • 

Henceforward, the notations of previous papers [HI HE] and the book [TT] 
are used. 

2 Integration of partial differential equations 
over Cayley-Dickson algebras. 

1.1. Notations and Definitions. Let A r denote the real Cayley-Dickson 
algebra with generators io, ■ «V-i such that io — 1, C- — —1 for each j > 1, 
ijik — —ikij f° r each j ^ fc > 1. It is supposed further, that a domain U in 
A™ has the property that 

(Dl) each projection Pj(U) =: Uj is (2 r — reconnected; 

(D2) ir s>Pt t(Uj) is simply connected in C for each k — 0, 1, 2 r_1 , s = 12k, 
P = *2fc+i) t G *4. r ,s, P an d u G C S)P , for which there exists z — u + 1 G Uj, 
where ej = (0, 0, 1, 0, 0) G A™ is the vector with 1 on the j-th place, 
Pj(z) = iz for each z G A™, z = jyjLi ^ ze ji ^ z e -A r for each j = 1, ...,m, 
m G N := {1,2,3,...}, 7r S)P)t (V) := {u : z G V, z = Eveb^vV, u = w s s + w p p} 
for a domain V in A r for each s ^ p G b, where t := Eveb\{s, P } ^v v G 
^r,s, P := {z G Ar : z = Eveb^vV, w s = w p = 0, w v G R Vv G b}, where 
b := {z , ^2^-1} is the family of standard generators of the Cayley- 
Dickson algebra A r - Frequently we take m = 1. Henceforward, we consider 
a domain U satisfying Conditions (Dl, D2) if any other is not outlined. 

The Cayley-Dickson algebra A r+ i is formed from the algebra A r with the 
help of the doubling procedure by generator i 2 r, in particular, Aq := R is 
the real field, A\ = C coincides with the field of complex numbers, A2 = H 
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is the skew field of quaternions, A3 is the algebra of octonions, A<± is the 
algebra of sedenions. The skew field of quaternions is associative, and the 
algebra of octonions is alternative. The multiplication of arbitrary octonions 
£,77 satisfies equations (1,2) below: 

(1) (£77)77 = £(7777), 

(2) £(£77) = mv , 

that forms the alternative system. The algebra A r is power associative, that 
is, z n+m = z n z m for each n, m G N and z G A r , it is non-associative and 
non-alternative for each r > 4. 

The Cayley-Dickson algebras are *-algebras, that is, there is a real-linear 
mapping A r 3 a 1— > a* G A r such that 

(2) a** = a, 

(3) (ab)* = b*a* for each a,b G A r . Then they are nicely normed, that is, 

(4) a + a* =: 2Re(a) G R and 

(5) aa* = a* a > for each 7^ a G A r . The norm in it is defined by the 
equation: 

(6) \a\ 2 := aa*. 

We also denote a* by a. Each non-zero Cayley-Dickson number 7^ a G 
^4 r has the multiplicative inverse given by ar 1 = a*/\a\ 2 . 

The doubling procedure is as follows. Each z G A r +i is written in the form 
z = a + 61, where l 2 = —1, 1 ^ ^4 r , a, 6 G ^4 r . The addition is componentwise. 
The conjugate of a Cayley-Dickson number z is prescribed by the formula: 

(7) z* := a* - 61. 

The multiplication is given by Equation 

(8) (a + /31)( 7 + 5\) = (a 7 - 5/3) + (5a + ^7)1 

for each a, /?, 7, <5 G A r , £ := ct + (31 G «4 r +i, 77 := 7 + 51 G »4 r+ i. 

The basis of A r+ \ over R is denoted by b r+1 := b := {1, z x , i 2 >-+ 1 -i}, 
where i 2 = — 1 for each 1 < s < 2 r+1 — 1, i 2 r := 1 is the additional element 
of the doubling procedure of A r +i from A r , choose ?2 r +m = im\ for each 
m — 1, 2 r — 1, i :— 1. 

1.2. Operators. An R linear space X which is also left and right A r 
module will be called an A r vector space. We present X as the direct sum 
X = X i © ... © X 2 r-ii2r-i, where X ,...,X 2 r_i are pairwise isomorphic real 
linear spaces. Certainly, for r = 2 this module is associative: (xa)b = x(ab) 
and (ab)x = a{bx) for all x G X and a, b G H, since the quaternion skew field 
A2 = H is associative. This module is alternative for r = 3: (xa)a = x(a 2 ) 
and (a 2 )x = a(ax) for all x G X and a G O, since the octonion algebra 
O = ^4 3 is alternative. 

Let X and Y be two R linear normed spaces which are also left and right 
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A r modules, where 1 < r, such that < ||ax||x = MIMIx = II^Hx and 
\\% + y\\x < \\ x \\x + \\y\\x f° r an x,y E X and a G A r . Such spaces X and F 
will be called A r normed spaces. 

We say that an A r vector space X is supplied with an A r valued scalar 
product, if 

{f,9) = T%£o(fj,9k)i* j ik, 
where / = foio + ...+ f 2^-112^-1, f,9 £ X, fj, gj G Xj, each Xj is a real linear 
space with a real valued scalar product, (Xj, (*, *)) is real linear isomorphic 
with (X fe , (*,*)) and (fj,g k ) G R for each j, fc. The scalar product induces 
the norm: ||/|| := ^fj)- 

An A r normed space or an A r vector space with A r scalar product com- 
plete relative to its norm will be called an A r Banach space or an A r Hilbert 
space respectively. 

We put X® k := X Cg) R ... ®r X to be the k times ordered tensor product 
over R of X. By L q ^ k (X® k , Y) we denote a family of all continuous k times 
R poly-linear and A r additive operators from X® fc into Y. If X and Y are 
normed A r spaces and Y is complete relative to its norm, then L qyk (X® k , Y) 
is also a normed R linear and left and right A r module complete relative to 
its norm. In particular, L q:1 (X,Y) is denoted also by L q (X,Y). 

If A G L q (X,Y) and A(xb) = (Ax)b or A(fcr) = b(Ax) for each x G X 
and 6 G A r , then an operator A we call right or left *4 r -linear respectively. 

An R linear space of left (or right) k times A r poly-linear operators is 
denoted by L lik (X® k ,Y) (or L r . :k (X® k ,Y) respectively). 

An R-linear operator A : X — > X will be called right or left strongly A r 
linear if 

A(xb) = (Ax)b or 

A(bx) = b(Ax) for each x G X and b G A r correspondingly. 

1.2.1. Examples. Each right A r linear operator A : X — > X such that 
Ax G X for each x G Xq is strongly right Ar- linear, since 

A(xb) = Y.J,k=QA(x j i j b k i k ) = Y.J, k =Q{Ax j )i j i k b k 

= Ef^JoiAxjiAikh = (Ax)b, 
since Axj G Xj up to an R linear isomorphism for each x = Xoio + ... + 
Xyr-iiyr-i ^ X and b = b i + ... + J^-i^-i G ^4 r with Xj G and bj G R 
for each j = 0, 2 r — 1. 

Particularly, if A G L r (X, Y) is right ^4 2 linear with X and F over A 2 , 
i.e. for X and F over the quaternion skew field H, then A is strongly right 
H linear, since the quaternion skew field is associative and x(ab) = (xa)b for 
each x G X and a, b G H. 

1.2.2. Lemma. If an invertible operator A is either right A r linear or 
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right strongly A r linear, then A 1 is such also. 
Proof. This follows from the equalities: 

(A-^b = (A- l Ax)b = xb = A-\A(xb)) = A-\(Ax)b) = A'^yb), 
where Ax = y, either x G Xq or x G X correspondingly and b G A r , since 
A(xb) = (Ax) b for each x G X Q oriel respectively and b G A r . 

1.3. First order partial differential operators. We consider an 
arbitrary first order partial differential operator a given by the formula 

(1) <rf = Y^ 1 *Wf/ d m)^ 

where / is a differentiable ^4 r -valued function on the domain U satisfying 
Conditions 1.1(D1, D2), 2 < r, i , i 2 7 '— i are the standard generators of 
the Cayley-Dickson algebra A r , ipj are real constants so that J2j ipj > 0, 
£ : {0, 1, 2 r — 1} — > {0, 1, 2 r — 1} is a surjective bijective mapping, i.e. 
£ belongs to the symmetric group (see also §2 in [16]). 

Making the substitution / \-t fi 2 r and using the embedding A r e — >■ A r+ \ 
it is always possible to relate the case of the operator a with ip ^ and a 
with tpj only for some j > 0, where is the doubling generator. 

For an ordered product {if---kf}q(k) °f differentiable functions s / we put 

(2) V{ 1 /... fe /} 9(fc) = E'] = oij{if-(dsf/dz^ j) )... k f} q(k) ip :j , 

where a vector indicates on an order of the multiplication in the curled 
brackets (see also §2 [HUE]), so that 

(3) Cr{lf--kf}q(k) = Y,s=l S °~{lf---kf}q{k)- 

Symmetrically other operators 

(4) tf^Zfjo^df/dz^i^, 

are defined. Therefore, these operators are related: 

(5) (afy = a(r). 

Operators a given by (1) are right A r linear. 

2. Integral operators. We consider integral operators of the form: 

/■oo 

(1) K(x,y) = F(x,y) + a F(z, y)N(ar, z, y)dz, 

J X 

where a is an R-linear partial differential operator as in §1.3 and a j is the 
non-commutative line integral (anti-derivative operator) over the Cayley- 
Dickson algebra A r from §§2.5 and 2.23 [16] or 4.2.5 [11] . where F and K are 
continuous functions with values in the Cayley-Dickson algebra A r or more 
generally in the real algebra Mat s (A r ) of s x s matrices with entries in A r . 
For definiteness we take a/dx J g(z)dz = J g(z)dz and a j g(z)dz calculated 
with the help of the left algorithm 2.6 [TO] 19] or 1.2.6(LJ) [JI]. In this case 
a J g(z)dz is the A r right linear operator in accordance with Lemma 1.2.2. 
If 7 : [a, b] — > A r is a function, then 
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(z) := supp | 7 (t i+ i) - 7(^)1 

is called the variation of 7 on the segment [a, b] C R, where the supremum 
is taken by all finite partitions P of the segment [a, b], P = {to = a < t\ < 
... < t n = b}, n G N. A continuous function 7 : [a, b] — > A r with the finite 
variation V a 6 7 < 00 is called a rectifiable path. 

Let a domain U be provided with a foliation by locally rectifiable paths 
{7° : a G A}, where A is a set (see below). We take for definiteness a 
canonical closed domain U in A r satisfying Conditions 1.1(D1, D2) so that 
00 G U, where A r = A r U {00} denotes the one-point compactification of A r , 
2 < r < 00. 

A path 7 :< a, 6 >— > U is called locally rectifiable, if it is rectifiable 
on each compact segment [c, e] C< a,b >, where < a,b >= [a,b] := {t G 
R : a < t < b} or < a, b >= [a, b) := {t G R : a < t < b} or < a, 6 >= 
(a, 6] := {t G R : a < t < b} or < a, 6 >= (a, 6) := {t G R : a <t < b}. 

A domain U is called foliated by locally rectifiable paths {7° : a G A} if 
7 :< a a , 6c >— )• U for each a and it satisfies the following three conditions: 

(Fl) UaeA7 a (< a m b a >) = U and 

(F2) 7 a (< a a , b a >) n 7 /3 (< ap, b p >) = ® for each a^eA. 
Moreover, if the boundary dU = cl(U) \ Int(U) of the domain U is non-void 
then 

(F3) dU = (UqgAi TM) U (U^a, 7 /3 (^)), 
where Ai = {a G A :< a a ,6/3 >= >}, A 2 = {a G A :< a a ,bp >=< 

a a ,bp]}. For the canonical closed subset U we have cl(U) = U = cl(Int(U)), 
where cl(U) denotes the closure of U in A v and Int(U) denotes the interior 
of U in A v . For convenience one can choose C l foliation, i.e. each 7° is of 
class C 1 . When U is with non-void boundary we choose a foliation family 
such that U a gA7( a a) = 9Ui, where a set dU\ is open in the boundary dU 
and so that w\qu 1 would be a sufficient initial condition to characterize a 
unique branch of an anti-derivative w(x) = T a f(x) = a f* x f(z)dz, 
where x G dUi, x G U, 7 a (to) — o x , l a if) — x f° r some a G A, t and 
t G< a a , b a >, 

^ J x f(z)dz = a J 7 «| [toit] /(z)dz. 

In accordance with Theorems 2.5 and 2.23 [16] or 4.2.5 and 4.2.23 [11] 
the equality 

(2) a x a ^ x g{z)dz = g(x) 
is satisfied for a continuous function gona domain U as in §1 and a foliation 
as above. 

Particularly in the class of A r holomorphic functions in the domain sat- 
isfying Conditions 1.1(D1, D2) this line integral depends only on initial and 
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final points due to the homotopy theorem [TD1 E] • 

We denote by P = P(U) the family of all locally rectifiable paths 7 :< 
a 7 , 6 7 >— > U supplied with the family of pseudo- metrics 

(3) p^( 7 , W ) := | 7 (a) -w(c)| + inf V a 6 ( 7 (t) - w(0(t)) 

where the infimum is taken by all diffeomorphisms : [a, 6] — » [c, d] so that 
0(a) = c and 0(6) = d, a < b, c < d, [a, b] C< a 7 , 6 7 >, [c, d] C< a w , 6 W >. 
We take a foliation such that A is a uniform space and the limit 

(4) ]im^ a p a ' b ' a ' b (^,j a ) = 
is zero for each [a, b] C< a a , b a >. 

For example, we can take A = R n_1 for a foliation of the entire Cayley- 
Dickson algebra A r , where n = 2 r , t G R, a G A, so that UaeA7 a ([0 ) °°)) = 
{z G Ar : Re{(z-y)v*) > 0} and U aG A 7 Q ((-oo, 0]) = {z£ A : i2e((* - 
2/)^*) < 0} are two real half-spaces, where v, y G A r are marked Cayley- 
Dickson numbers and t) / 0. Particularly, we can choose the foliation such 
that 7 a (0) = y + a±Vi + ... + a^-i^V-i and 7°(t) = £i>q + 7 a (0) for each 
t G R, where t>o, ••.,"U2 r -i are R-linearly independent vectors in A r . 

Therefore, the expression 

/•oo r 

(5) a I g(z)dz ■= a g(z)dz 

denotes a non-commutative line integral over A r along a path 7 so that 
l a {tx) — % and lim{,_>6 a 7 a (t) = 00 for an integrable function g, where t x G< 
a a ,6 Q >, a G A, a Q = a 7 <*, 6 Q = & 7 «. It is sometimes convenient to use the 
line integral 

(6) fj / g(z)dz •■= * g(z)dz, 

Jx h a \< aa ,t x] 

when lim a ^ aa ^ a (t) = 00. 

If a J g(z)dz = u~ l g corresponds to the right (or left) algorithm of the 
non-commutative line integration over the Cayley-Dickson algebra A r , then 
H T [{a I g( z )dz)i2r] corresponds to the left (or right correspondingly) algo- 
rithm such that 

M(a/)<5r] = Y%Jo l {dr/dz m )i^ = *(/*), 
where / is a differentiable A r valued function. 

For example, one can take a function g(z) as F(z, y)K(x, z) depending on 
additional parameters x,y,t, .... 

A function F has the decomposition 

(7) f = e?: 1 f^-, 

where each function Fj is real-valued for a function F. 

Put for convenience a = I, where I denotes the unit operator, a m 
denotes the m-th power of a for each non- negative integer < m G Z. 
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3. Proposition. A family T> r of all differential operators with constant 
A r coefficients is a power associative real algebra with a center Z(T> r ) consist- 
ing of all differential operators with real coefficients and with a unit element 
I. 

Proof. Let A be a differential operator with constant A r coefficients, 
then it can be written in the form: 

(1) Af = j:f-J [A jis7 r s f]i*, 

where Aj >8 is a differential operator with real coefficients for each j: 

(2) Aj jS ^2ct ^j,s,ad , 

aj tSt a G R for each a = (a , a^-i), d a = aQ a ' a '^-i > l a l = a o + ... + 
a^-i, ak = 0,1,2,... is a non-negative integer for each k = 0,...,2 r — 1, 
iT s : X — > X s denotes the R-linear projection operator (see §1.2), 

(3) Af = vj: „(.t,,./;)/; 

for each / in a domain Dom(A) of the operator A, where 

(4) / = ^ /A, 

f s is a real-valued function for each s, f s = n s f. If P n {z) is a polynomial 
with A r coefficients of the variable z G A r (see also §2.1 0, [ID]), then 

(5) AP n ( Z ) = Zf-Jo^aM<n a i,sAd a (Pn)s(z)}i*, 

consequently, P n belongs to the domain Dom(A) of the operator A. Thus 
Dom(A) is non-void for each differential operator A G T> r . From Formulas 
(1 — 4) we get, that the sum A + B and the product BA of any differential 
operators A,B eV r belongs to T> r , since 

(6) BAf = B(Af) = ESio(-l) si9n()) (5y4^)^ = Cf with 

(7) C ktS = ^(-ly^B^, 

where sign(x) = for x — 0, sign(x) = 1 for x > 0, sign(x) = —1 for 
x < 0. Certainly L4 G P r for each t G R and A G P r . From Formulas 
4(7 - 9) [T6] or 4.2.4(7 - 9) [H] we deduce that 

(8) A 2 / = A(Af) = Y$£{A)f)ij + ^U'UoAj + AjA )f]i* 

for each / G Dom(A 2 ), since Jx = X^o* ^ji^h f° r each i6l, where 

(9) .1,/ X: ,/ ^y,j>- 

since (i*) 2 = i 2 . Therefore, by induction we infer that A 2k = (A 2 ) h and 
A 2k+1 = (A 2 ) k A = A(A 2 ) k for each natural number k, where A = I denotes 
the unit operator, i.e. If = f for each function / : A r — > A r . Thus the 
algebra V r is power associative: 

(10) A k A m = A k+m for each natural numbers k, m G N and every A G T> r . 
Moreover it contains the unit element I. 

If A G Z(T> r ), i.e. an operator A commutes with every differential oper- 
ator B G D r , [A,B] := AB — BA = 0, then it commutes with all generators 
{ij, d Xk i* : j, k = 0, 2 r — 1} of the algebra V r . From Formulas (1 — 4, 6, 7) 
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it follows that %, s , a = for j > 0, where %, s , a £ R for each j, s and a. That 
is, all the coefficients aj :Sja i* of the operator A are real, since i — 1. 

3.1. Corollary. Lei A G D r k a differentiable operator with constant 
A T coefficients and n : X — >■ X fre £/ie R linear projection (see §1.2), then 
[(I — ir )A] 2k is with real coefficients for each natural number k e N. 

Proof. In view of Formulas 3(1, 9) the equality is valid: 

(1) (/ - n Q )Af = Y.fJiiAjf)^, consequently, 

(2) [{I-^)AYf = YZ = \\A)f)il 

That is, [(/ — 7r )A] 2 is the differential operator with real coefficients, since 
ij G R for each j. 

3.2. Remark. Mention that in the considered case vp = and the 
operator a 2 is with real coefficients due to Corollary 3.1, hence the square of 
the ant i- derivative operator ( CT J x ) 2 or a~ 2 is also with real coefficients, when 
each ipj is a real constant and ip = 0. 

4. Proposition. Let 

(1) Jim V 2 fc2 < 2 a:F(z,y)K(x,z)=0 

for each x,y in a domain U satisfying Conditions 1.1(D1, D2) with oo e U 
and every non-negative integers < k, s, n e Z suc/i that k + s + n < m. Sup- 
pose also that a f£° d%d!j}d£[F(z, y)K(x, z)]dz converges uniformly by parame- 
ters x,y on each compact subset W C U C A 2 for each \a\ + \[5\ + \oj\ < m, 
where a = (a , o^-i), \a\ = a + ... + ct2 r -i, d x = /8xq° ...dx^J] 1 . 
Then the non- commutative line integral c f£° F(z,y)K(x, z)dz from §2 satis- 
fies the identities: 

roc roc 

(2) o™ a \ F(z,y)K(x,z)dz = 2 o™ c / F(z, y)K(x, z)dz+A m (F, K)(x, y), 

Jx Jx 

POO POO 

(3) V 2 m CT / F(z,y)K{x,z)dz = {-l) m2 a? <T F(z,y)K(x, z)dz+B m (F,K)(x,y), 

Jx Jx 

where 

(4) A m (F,K)(x,y) = - 2 a™~ 1 [F(x,y)K(x, z)}\ z=x + a x A m _i(F, K)(x,y) 
for m > 2, 

(5) B m (F,K)(x,y) = (-l) m 2 a™~ 1 F(x,y)K(x, z)\ z=x +[ 1 a z B m _ 1 (F(z,y), K(x,z))]\ z=; 
for m > 2, 

(6) A 1 (F,K)(x,y) = B 1 (F,K)(x,y) = -F(x,y)K(x,x), 
a x is an operator a acting by the variable x G U C A r . 

Proof. From the conditions of this proposition and the theorem about 
differentiability of improper integrals by parameters we have the equality 
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* r d£dl% [F(z, y)K(x, z)]dz = ^ (7 £° dt [F(z, y) K(x, z)]dz 
for each |a| + + |a>| < m (see also Part IV, Chapter 2, §4 in [5J). 

In view of Theorems 2.5 and 2.23 and Corollary 2.6 [16] or 4.2.5 and 
4.2.23 and Corollary 4.2.6 [1JJ the equalities 



poo 

(7) o x a \ g{z)dz = -g(x) and 

Jx 



X 



(8) (j / [a z f(z)]dz = f(x) - f( x) 

are satisfied for each continuous function g and continuously differentiable 
functions g and /, where is a marked point in U, 

pOO z 1 pOO 

(9) 1 a zrT / F(2,y)K(a;, z)c?2: := ^ <r / Oj[(^F(^, y)/dz^ j) )K{x, z)]tpj}dz and 

2 r -l 



px " * px 

(10) 2 o za F(z,y)K(x,z)dz := J2 ° {i*j[F(z, y)(dK(x, z)/dz^ j - ) )]ip j }dz &nd 

J OX j = Q J OX 
pOO z 1 pOO 

(11) 2 a xa I F(z,y)K(x,z)dz := J2 {i*j[^(z,y)(dK(x, z)/dx^ j) )]^ j }dz. 

JX • q JX 

Therefore, from Equalities (7,8), 1.3(3) and 2(5) and Condition (1) one gets: 

poo poo 

(12) a xa F(z,y)K(x,z)dz = 2 a x a F(z, y)K(x, z)dz- F(x, y)K(x, x), 

Jx Jx 

since F(z,y)K(x,z)\™ = — F(x, y)K(x, x), that demonstrates Formula (2) for 
m = 1 and A\ = —F(x, y)K(x, x). The induction by p — 2, m gives: 

poo 

(13) alJ F(z,y)K(x,z)dz = 

Jx 

poo 

a x ?v p x~ 1 ° \ F{z,y)K{x,z)dz] + a x A p _i(F, K){x, y) 

Jx 

poo 

= 2 a p xa F(z,y)K(x,z)dz 

Jx 

+ 2 a p x - 1 A 1 (F(x,y),K(x,z)\ z=x ) + a x A P ^(F, K)(x,y). 

Substituting in these expressions Ai one gets Formulas (2, 4) for each m > 2. 
Then from Formulas (7,8) and Condition (1) we infer also that 



poo poo 

(14) 1 a za F(z,y)K(x,z)dz = - 2 o z / F(z, y)K(x, z)dz+F(z, y)K(x, z) 

Jx Jx 

poo 

= -F(x,y)K(x,x) - 2 a z a / F(z, y)K(x, z)dz, 

Jx 

that gives Formulas (3) for m = 1 and (6). Then we deduce Formulas (3,5) 
by induction on p = 2, m: 

poo 

(15) l a p z<7 F(z,y)K(x,z)dz = 

Jx 
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(-l)^ 1 V 2 [Vr 1 (T / F(z,j/)K(a;,*)tfe] + V, B p _i(F(z,y),KOr,z))| z=a; 

J X 

roo 

= {-iy 2 al„ F(z,y)K(x,z)dz 

J X 

+ 2 ar 1 F(z,y)K(x,z)\ z=x + 1 a z B p ^(F(z, y), K(x,z))\ z = x . 

4.1. Remark. The center of the Cayley-Dickson algebra A r for r > 2 is 
the real field R. 

Let Mat s (A r ) denote the left and right A r module of square sxs matrices 
with entries in A r . It is possible to consider the quaternion skew field Hj^,l 
with generators J, K, L realized as 4 x 4 square real matrices. Putting ijY = 
{yk,flj) — Yij for each real matrix Y with elements y^i G R for each k, I 
and every generator ij of the Cayley-Dickson algebra we naturally obtain 
the quaternionified algebra (A r )n JKL (see also |8j). On the other hand, 
H C Mat 4 (R), consequently, Mat s ((A r )Hj K L ) C Mat is (A r )- 

5. Corollary. If suppositions of Proposition 4 are satisfied, then 

(1) A 2 {F,K)(x,y) = -a x [F(x,y)K(x,x)]- 2 a x [F(x, y)K{x, z)]\ z=x , 

(2) A 3 (F, K)(x,y) = -a 2 x [F(x,y)K(x,x)\ 

-Vx( 2 o- x [F(x,y)K(x,z)]\ z=x ) - 2 a 2 x [F(x,y)K(x, z)]\ z=x , 

(3) B 2 (F, K)(x, y) = - 1 a x [F(x, y)K(x, x)] + 2 a z [F(x, y)K(x, z)]\ z=x , 

(4) B 3 (F,K)(x,y) = - l cr 2 x [F(x,y)K(x,x)] 
+ 1 a x ( 2 a z [F(x,y)K(x,z)]\ z=x ]) - 2 a 2 [F(x, y)K(x, z))\ z=x . 

(5) A 2 (F, K)(x, y) - B 2 (F, K)(x, y) = -2 2 a x [F(x,y)K(x,x)\, 
where a x K(x, x) = [a x K(x, z) + a z K(x, z)]\ z=x , 

(6) A 3 (F, K){x, y)-B 3 {F, K)(x, y) = -(3 2 a 2 x + 2 a x 2 a z +2 2 a z 2 a x )[F(x, y)K{x, z)} 
-(2 l a x 2 a x + 2 a x V x )[F(a:,y)K(x,x)]. 

Particularly, if either p is even and ipo = 0, or F G Mat s (R,) and K G 
Mat s (A r ) , then 

(7) 2 aP[F(z,y)K(x,z)} = F(z,y)<r*K(x,z) and 2 a? [F(z, y)K(x,z)] = F(z,y)aPK(x 
Proof. Formulas (1 — 6) follow from Equalities 4(4 — 6). In particular, if 

p is even and ipo = 0, p = 2k, k G N, then in accordance with Corollary 3.1 
and Formulas 4(7 - 9) [IS] or 4.2.4(7 - 9) [TTj 

a*f{x) = A k f(x) 
for p times differentiable function / : U — > A r , where 

Af = Ei bjd 2 f(x)/dx 2 , b 3 = i 2 . 1(j) e R. 

Evidently the operators 2 a v x and 2 o~ v z commute with the left multiplication 
on F(z,y) G Mat s (R) in accordance with §4.1, that is 2 a x '[F(z, y)K(x, z)] = 
F(z, y)<r>K(x, z) and 2 al [F(z, y)K(x, z)\ = F(z, y)a*K(x, z). 

6. Remark. Expressions of functions A m and B m depend not only on F 
and K, but also on a, that is on coefficients ipj for j = 0, 2 r — 1 with 2 < r 
in accordance with 554 and 5. 
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7. Method of non-commutative integration of vector partial 
differential equations. 

Let an equation over the Cayley-Dickson algebra A r be given in the non- 
commutative line integral form: 



where K, F and N are continuous integrable functions as in §2, p e R \ 
{0} is a non-zero real constant. These functions may depend on additional 
parameters t,r, .... 

The first step consists in a concretization of a function N and its expression 
throughout F and K. Suppose that an operator 

(2) (I — A x )F(x : y) = K(x, y) is invertible, so that (I — A^,) -1 is continuous, 
where / denotes the unit operator, while 



is an operator acting by variables x. 

On the second step two R-linear differential or partial differential opera- 
tors L s over the Cayley-Dickson algebra A r are given: 



where / is a differentiable function in the domain of L s , L s jg is real- valued 
function for each ord(L s ) times differentiable real-valued function g in the 
domain of L s j for every j (see also Formulas 3(1 — 4)), where ord(L s ) de- 
notes the order of the partial differential operator L s . On a function F the 
conditions either: 



for each s = 1, 2 and 1 < k < m are imposed, where jC are constants jC G A r , 



^ k C {0, 1, T - 1} for each k, U fc ^ = {0, 1, 2 r - 1}, V k n % = for 



each k y£ I, 1 < m < 2 r . Coefficients jC or operators L s j may be zero for 
some j. 

On the third step a function K is calculated from Equation (2). 
Acting by the operator L s from the left on (2) and using Conditions either 
(4) or (5) one gets either 

(6) L S [(I — A X )K] = or the equalities: 

(7) E je * k ij{ jcL a ,o[(/-A !e )K]+L aJ [(/-A x )K]} = for each k = l,..,m. 
Therefore, using Conditions either (4) or (5) we infer that 




P a iT F(z,y)N(x,z,y)dz 



(3) /-,/ v ';(/->,,/!. 
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(9) R s (f) = (I - A X )(LJ) - L S [(I - A x )f] 

is formed with the help of commutators [A, B] = AB—BA and anti-commutators 
{A,B} = AB + BA of operators (I - (A x ) ) , (A*)*, L 9j , fc,j = 0, 2 r -l. 
A function N and operators L\ and L 2 are chosen such that 

(10) R S (K) = (I - A x )M a (K) for s = 1, 2, 

where M S (K) are functionals of K which generally may be non-R-linear. In 
view of Condition (2) the function K must satisfy partial differential equations 

(11) L S K — M S (K) = for s = 1, 2, 

which generally may be non-R-linear. Thus each solution K of the R-linear 
integral equation (1) is also the solution of partial differential equations (11). 
Frequently particular cases are considered, when Equations (4) correspond to 
an eigenvalue problem for s = 1 and to an evolution in time problem for s = 2. 
Generally operators L sQ , L Si 2 r -i can be chosen A r vector independent and 
jC G Rifc(j) for each j and s = 1, 2 with fc = G {0, 1, .., 2 r — 1}. 

The Euclidean space R 2r is the real shadow of the Cayley-Dickson algebra 
A r , that is, by the definition A r considered as the R linear space is isomorphic 
with R 2 \ The Lebesgue (non-negative) measure /i on the Borel a-algebra 
B(R 2r ) of the Euclidean space R 2r induces the Lebesgue measure on B(A r ). 
Therefore, the Hilbert space X = L 2 (A r , /i, A r ) of all Lebesgue measurable 
functions / : A r — > A r with integrable square module \f\ 2 , i.e. 

{HN) || /|| 2 := J Ar \f(z)\ 2 fi(dz) < oo, and with the A r valued scalar 
product 

(SP) (f,g):=f Ar f*(z)g(zMdz) 
exists. Analogously the Hilbert space X = L 2 (A r , n, Mat s (A r )) is defined 
with J2j,k=o fj,k( z )9j,k(z) instead of f*(z)g(z) in the integral in Formula (SP), 
where fj^ G A r denotes a matrix element at the intersection of row j with 
column k of a matrix / G Mat s (A r ). 

Let A be an R linear A r additive operator A : D(X) — > Y, where D(A) is 
a domain of A dense in X, D(A) C X, X and Y are Hilbert spaces over A r . 
Then its adjoint operator A* is defined on a domain consisting of all those 
vectors y EY such that for some vector z G X the equality (x, z) = (Ax, y) 
is valid for all x G D(A). For such y G Y we put A*y = z. If A* = A we 
say that A is self-adjoint. Thus (Ax,y) = (x,A*y) for all x G D(A) and 

y e D(A*). 

If an operator A is self-adjoint A = A*, then (Ac,y) = (x,Ay) for all 
x,y G D(A). 

When D(A) is dense in X and (Ax, y) = (x,Ay) for all x, y G D(A), we 
say that A is symmetric. A self-adjoint operator is maximal symmetric. 
We define a self-adjoint operator A to be positive, when (Ax, x) > for 
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each x G D(A). 

An R linear A r additive operator A is called invertible if it is densely 
defined and one-to-one and has dense range 11(A). 

The operator A* A is self-adjoint and positive. If A is invertible, then 
(A* A)- 1 A* C A- 1 (see also [31 [T5]). 

If an expression of the form 

(12) £ fc [(I - A x ) k f(x, y)\ k g(y) = u(x, y) 

will appear on a domain U, which need to be inverted we consider the case 
when 

(RS) (I — A x ) is either right strongly A r linear, or right A r linear and 
kf G X for each k, or R linear and tg(y) G R for each k and every y E U, 
at each point a? G U. 

When an operator (I — A x ) is invertible and Condition (RS) is satisfied, 
Equation (12) can be resolved: 

(13) J2k kf(x, y) k g(y) = (I - A x )~ l u(x, y). 

If condition (RS) is not fulfilled, the corresponding system of equations 
in real components (A x )j >s , kfs and kg s can be considered. 

If A : X — > X is a bounded R linear operator on a Banach space X with 
the norm ||A|| < 1, then 

The anti-derivative operator g \-t a f® x g(z)dz is compact from C°(V, A r ) into 
C°(V,A r ) for a compact domain in A r , where C°(V,A r ) is the Banach 
space over A r of all continuous functions g : V — > A r supplied with the 
supremum norm \\g\\ := swp xeV \g(x)\, ox is a marked point in V, x G 
V. Therefore, an operator A x will be invertible with a suitable choice of a 
function F satisfying the system of R linear partial differential equations (4) 
or (5) and a real non-zero parameter p ^ 0. 

7.1. Proposition. Le£ V be a compact domain in the Cayley- Dicks on 
algebra A r with 2 < r and let its foliation be with rectifiable paths and satisfy 
Condition 2(4) and A be a compact subset in R 5 ! . Then the anti-derivative 
operator a f x from §2 is compact from C°(V, A r ) into C°(V, A r ). 

Proof. As usually C l (V,A r ) denotes the space of all continuously differ- 
entiable functions / : V — > A r with the supremum norm 

\\f\\ c i := sup xeV \f(x)\ + E?Lo lsu Pxey \df(x)/dx j \. 

The decomposition is valid: C S (V, A r ) = C S (V, R)z © ... © C S (V, R)i 2 - -i 
for s = 0, 1, .... On the other hand, the partial differential operators C 1 (V, R) 3 
f I—)- i^(df(x)/dxj) G C°(V, R) are Rlinearly independent, since Re(iji* k ) = 
for each j ^ k. In view of Theorems 2.6 [9j [10] and 1.2.7, 4.2.5 and 
4.2.23 and Corollary 4.2.6 [UJ the ant i- derivative mapping C°(V, A r ) 3 g i— > 
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a J x g(z)dz G C 1 (V, A r ) is continuous, since a x a j x g(z)dz = g(x). But the 
embedding C 1 (l / , A r ) C°(V,Ar) is the R linear Ar additive compact op- 
erator. Therefore, the anti-derivative operator a J x : C°(V,A r ) — > C°(V, A r ) 
is compact. 

8. Example. Let 

(1) N(x,z,y) = K(x,z) in 7(1). 

Acting from the left on both sides of Equation 7(1) one gets: 

(2) L s K(*,y) = L s pair f(z,y)K(x,z)dz. 
Take the hyperbolic partial differential operator 

(3) L x = ,al - 2 *l 

of the second order, where the operator fc<r is with coefficients ki>j £ R f° r 
each j = 0, T — 1 and a transposition £ k G S 2 r (see §1.4). Condition 7(4) 
for s = 1 means that L X F = 0, i.e. 

(4) 1 a 2 z F(z,y) = 2 a 2 y F(z,y). 

Then due to Proposition 4 and Corollary 5 we deduce that 
(i<7* - 2 ^)K(x,y) = p(i(j2 - 2<T J) lff / x °°F(3,j/)K(a;,2)cfc 
= P K * L 00 F(z, y)K(x, z)d*+p iA 2 (F, K)(x, y)-p Jof 1<r XT F(z, y)K(x, 
= P(?^- i* £° F (^, y)K(x, z)d*+p 1 A 2 (F, K)(x, y)-p 1 5 2 (F, K)(x, y), 

consequently, 

(5) (/- A^td^ - 2S 2 )K(x,y)] = p !A 2 (F,K)(x,y) - p i-B 2 (F, K)(x,y) 
= -2p ?<7 x [F(a;,j/)K(a;,a;)] = -2p ^{[(/-As)K(z,y)]K(rr,:r)}, 

where the terms fcA s and kB s correspond to ^a. 

If the Cayley-Dickson algebra A r is either with r = 2 and F G Mat s (R) and 
K G Mat s (A), s G N, or F G R and Kei 3 with r = 3, then 

? < T a: {[(/-A a: )K(x,y)]K(x,a;)} = [(/ - A x )K(z, y)][ ^(a;, z)], 
since R is the center of the Cayley-Dickson algebra A r with r > 2, < 
F(z, y), K(x,z), i<j x K(x,x) >= in these cases, where < a,b,c >:= (ab)c — 
a(bc) denotes the associator of Cayley-Dickson numbers a,b,c G A r . Then 
Conditions 7(1, 2) imply that the function K satisfies the non-linear partial 
differential equation: 

(6) (kt 2 x - 2 <j 2 y )K(x,y) + 2pK(x,y)[ 1 <j x K(x,x)}=0. 

If put u(x) = 2 icr x K(x,x) over the quaternion skew field H = A 2 , i.e. 
for r = 2, and substitute K(x, y) = k) exp(JRe(ky)) into (6), we deduce 
that a function <3> satisfies Schrodinger's equation: 

(7) 1 a^(x, k) + fc)(p« + Ej kjtj 2 tf) = 0, 

where k G H, since H is associative and the generator J commutes with 
io, i 2 r-i, also exp(JRe(ky)) commutes with $(x,k). 

Now we take the third order partial differential operator with iip = 2 ipo = 

0: 
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(8) L 2 f = ( 3 a t +i^+3 2 a y x o 2 x + 3 2 a 2 i<j x + 2 <rJ)/- Then we put 

(9) L 2j 7 = [ 3 ^ 3 tfA + d^j 9 x ei0 )+ 3 rtAtoV)) i f7 x+(2'0j ^ 2(3) + 3 i^A €l0) ) 2^]/ 
and impose Condition 7(5) for s = 2: 

(10) L 2 jF(x,y) = for each j, 

where the operator 3 <7 t is with real constant coefficients 3 ipj and a trans- 
position £ 3 G SV- We suppose that functions F and K may depend on t. 
Therefore, we get from Equation (2) and Condition (10): 

(11) ( 3 <7i + x al + 3 2 a y \o 2 x + 3 2 a 2 x a x + 2 a*)K(x,y) 

= p( 3 a t +i(T^+3 2 a y x a 2 x + 3 2 a 2 y x a x + 2 aJ) 1<T /~ F(z, y)K(x, z)cfe = 

-p(K 3 + 3 fo 2 + 3 K 3 ) i^F(^)^,^+ 

P(i<7t + i(r x + 3 2 o y x a 2 x + 3 2 <j 2 y x a x + 2 <rJ) 1( , J™ F(z,y)K(x, z)dz = 
I + pfrt 1 aI x x F(z,y)K(x,z)dz, 
where I — h + I 2 + I 3 , 

(12) h = p( l0 f- ^ 3 ) jrF( Z)!; )K(v)i 

= PM+ ?^ 3 ) 1CT r F(z, y)K(x, z)dz+p iA 3 (F, K)(x, y)-p 1 B 3 (F, K)(x, y) 

(13) I 2 = 3p( 2 a y x a\ - 2 a y \a 2 z ) ia F(z, y)K(x, z)dz 

= 3p 2<T y [(l<T 2 x - \a 2 z ) ia C F(z, y)K(x, z)dz+ ± A 2 (F, K)(x, y)- 1 B 2 (F, K)(x, y)}, 

(14) I 3 = 3p( 2 a 2 x a x - 2 a 2 y \a z ) 1<r f~ F{z, y)K(x, z)dz 
= 3p 2° 2 y { \o x + \a z ) 1CT JT F(z,y)K(x,z)dz, 

since iAl = 1 Bi. Equations (4, 14), 4(3) and 5(1,3) imply that 

(15) I 3 = 3p \a 2 z ( \a x + \a z ) ia f™ F(z, y)K(x, z)dz 

= 3p( 2 a z -la x ){( 2 a x + 2 a z )[F(x,y)K(x,z)]}l^ 1<r / x °° F(z, y)K(x, z)dz. 

Particularly, if the Cayley-Dickson algebra A r is with 2 < r < 3 and 
F G Mat s (R) and K G Mat s (A), s G N for r = 2, s = 1 for r = 3, 
then Equations (6, 13) imply that 

(16) I 2 = -3 2 a y [K(x, y)-F(x, y)]u(x)+ 3p 2 (r y [iA 2 (F, K){x, y)- 1 B 2 (F, K)(x, y)}, 
since 

(17) -2p 1CT r F(z, y) [K(x, y)da x K(x : x))]dz = -p ia F(z, y) [K(x, z)u(x)]dz 
= ~P[ ia 1^ F(z,y)K(x, z)dz]u(x) = [K(x,y) - F(x,y)]u(x), 

where u(x) = 2 ia x K(x, x). We deduce from Formulas (6, 10, 11, 15 — 17) that 

(18) ( 3 a t + 1 a x +3 2 a yl a 2 x +3 2 al 1 a x + 2 a^)K(x,y)+3 1 2 a y [K(x,y)u(x)} = 

P(!*t + \°l + 3 fc* fc 2 + 3 \°l i°x + i<r/*°° F(z,y)K(x,z)dz 
+3p 1<T F(z, y)[la z K(x, z)u(x)]dz + T, where 

T = p ^(F, K)(x, y)-p 1 5 3 (F, K)(x, y)+ 3p 2 <T y [iA 2 (F, K)(x, y)- ± B 2 (F, K)(x, y)}+ 

3p la y [F(x, y)u{x)]+3p ( \o z 2 a x + \a 2 z - \o x \o x -\o x 2 a z )[F(x, y)K(x, z)]\ z=x . 

Then for each two continuously differentiable A r valued functions G(z,y) 
and K(x,z) one has 

(V,- 2 a x - 2 <j z )[G(z,y)K(x,z)] = ( 1 <j z + 2 a x + 2 a z )[G(z, y)K(x, z)] = 
(cr x + <T z )[G(z,y)K(x,z)], 
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where K(x,z) := K(— x, —z) for each —x, —z G U. Therefore, the identity 

[aLA 1 ^, 2 a x ]]{G(z,y)K(x,z)}dz = 
is satisfied, consequently, in the considered case 2 < r < 3 and F G Mat s (H) 
and K G Mat s (A r ), s G N for r = 2, s = 1 for r = 3, we get 

(19) [(/ - A x ), [V,, Vj]{F(x, y)K(x, x)} = 0, 

since R is the center of the Cayley-Dickson algebra A r and 

(20) [ 2 a x , 3 <j x ]F(z,y){F(x,z)K(x,x)} = F(z, 2 /)([ 1 ^, V x ]{F(x, z)K(x, *)}) 
due to Formulas (1), 4(7,8) and 7(2), since [V x + 2 a x , 1 a x — 2 a x \ = 
-2^, 2 a x ]{G(x,y)K(x,x)}. 

Substituting the expressions of A 2 — B 2 and A 5 — B 3 from Corollary 5 
and using Formulas (2, 6, 19, 20) one gets: 

(21) T = -p(3 2 a 2 x + la x 2 a z + 2 2 a z la x )[F(x,y)K(x,z)]\ z=x 
-p(2 \o x \o x + \a x \a x )[F(x,y)K(x,x)} 

+3(l-p) 2 (?y[F(x, y)u{x)]+3p ( \a z \a x + \a 2 z - \a x \a x -\o x \a z ) [F(x, y)K(x, z)]\ z=x 
= -3p \a x [F(x, y)u(x)]-3p(ja 2 £ - ja 2 )[F(x, y)K(x, z)]\ z=x +3(l-p) 2 a y [F(x, y)u(x)] 
+p[fc 2 , lvx}[F(x,y)K(x,z)}\ z=x + p[\a x , la x ][F(x,y)K(x,x)] 
= -3p \a x [F(x,y)u(x)] + 3pF(x,y)[K(x,x)u(x)}+3(l-p) 2 cr y [F{x, y)u(x)\ 
+p[fc 2 , \cr x ][F(x,y)\i(x,z)\\ z=x + p\ L 1 a x , \a x ][F(x,y)K(x,x)\ 

= -3p \a x [K(x, y)u{x)]+3p 2 i(J x {{ ia F(z, y)K(x, z)dz)u(r])}\ v=x +3pF(x, y)[K(x, x)u(x)} 
+3(l-p) 2 a y [F(x,y)u(x)]+p[ 2 a z , la x ][F(x,y)K(x, z)]\ z=x +p[\ <2xi i&x 

][F(x,y)K(x,x)] 

= -3p \a x [F(x,y)u(x)\ + 3p 2 2 a x [( lfT f x °° [(F(z, y)K(x, z))u{x)]dz 
+3(l-p) 2 a„[(/-A x )K(a;, J /)«(a;)]+p[(/-A e )K(a;,y)]{[i<T z , i<J x ]K(x, z)]\ z=x }+ 
p(J - A x ){[\a x , ja x ][K(x,y)K{x,x)]}, 

since 2 < r < 3, F G Mat s (R) and K G Mat s (A r ) with s G N for 
r = 2 and s = 1 for r = 3, hence < F(x,y),K(x,x),u(x) >= 0, where 
< a,b,c >:= (ab)c — a(6c) denotes the associator for each Cayley-Dickson 
numbers a, b, c G A r , since i ^fc — ^fc^o — 0, ^4 = — ikij f° r eac h j ^ k > 1, 

(22) [<r 2 , cr x ] = ((T 2 (T X - cr x a z ) = —[a x , a z \. 

We now take p = 1. Therefore, in accordance with Formulas (19 — 21) and 
7(13) the equality 

(23) { 3 a t +i^+3 2 a y l( r 2 x + 3 2 a 2 y x a x + 2 a*)K(x,y) 
+6(\a x + l<jy)[K(x,y)( 1 cr x K(x,x))]-K(x,y){[ 1 cr z , i<r x ]K(x, z)]\ z=x } 
-[Wx, i<r x ][K(a:,y)K(a:,a;)] = 

follows, when the operator (/ — A x ) is invertible. Particularly, for s = 1 and 
\<j = 2 o with = and 3 a t = d/dt differentiating Equation (23) with 
the operator \o x and then restricting on the diagonal x = y and taking into 
account Formulas (1,6), 7(1) and (19) one gets the equation 

(24) u t (t, x) + 6 ia x [u(t, x)u(t, x)} + i<j%u(t, x) = 

of Korteweg-de-Vries' type, where u(t,x) = 2 i<j x K(x, x), since a(f(x)g(x)) = 
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(V + *a)(f(x)g(x)), [V, 2 a](f(x)g(x)) = -±[<t, V - 2 a](f(x)g(x)), 
<x[V, 2 a](f(x)g(x)) = -[V, 2 (7]<t(/(x)^(x)) and hence a^a, 2 a](f(x)g(x)) = 
— [V, 2 a]a~ 1 (f(x)g(x)) for a = \O x with 1^0 = and continuously differen- 
tiable functions f(x) and g(x) (see Remark 3.2). 

8.1. Theorem. A solution of partial differential Equation (23) with 
iV'o — 2^0 — over £/ie Cayley-Dickson algebra A r with 2 < r < 3 is given by 
Formulas (2 — 4, 9, 10) with p = 1 whenever the appearing integrals uniformly 
converge by parameters on compact sub-domains (see Proposition 4) and the 
operator (I — A x ) is invertible and F G Mat s (R) and K G Mat s (.4 r ), s G N 
w/ien r = 2, s = 1 w/ien r = 3. 

9. Example. Consider the integral equation 

(1) K(x,y) = F(x,y) + f F(«,y)[F(z,«)K(x,«)](fa)d«, 
where p is a non-zero real constant. We take, for example, 

(2) L x = o x - x o y with = iV'o = 0. 
A solution of the equation 

(3) L 1 F(x,y) = 

has the form F(x, y) = G( (a '" >+ 2 (lQ ' i/> ) or we shall write F( (a ' 3!> + (ia ' y> ) instead 
of F(x, y), where 

(4) (a, a; >:= Vj „ ; ",-'•,/, 

for a,x <E A r , x = Y?j=o x jij, x j £ R- f° r each j, = for ipj ^ 
and a,j = 1 for ^ = 0. Using suitable change of real variables Xj,yj we 
can suppose without loss of generality, that ipj, \ipj G {0,1} for each j = 
0, T — 1. Therefore, a solution of Equation (3) can be written as 

(3.1) F( (M + y> ) as well, where bj = aj ^ G {0,1} for each j. Differ- 
entiable functions of the form F( ^ b ' x ^ v> ) are also solutions of the system of 
partial differential equations 

(3.2) L ltj F(&^) = 0, where 

(3.3) /.,..; = i>jd Xj - iiJjd yr 

It is convenient to introduce the notation 

(5) K 2 (x, z) := a J °° F(te^^)K(x, x + CR- 
Using Condition (3) we rewrite Equation (1) in the form: 

(6) K(x,y) = F( ^>+^y> ) 

+! *f5° o Jo°° [F( ( a '*+*>+b a <v> ) (F( (°.*+<;>+(i",*+v> )K(x,x + C))dC]dri 
= F ( (^>+(i^> ) + £ CT / °° F{ {a > x+v> + {ia > y> )K 2 {x,x + V )d V . 
Then let us put: 

(7) AJ(y) := f CT / °° a / °° [F( ^>+^v> ) (F( (°*+<>+i™+» )f( x +())d(]dri 
for a continuous function /, consequently, 

(8) (I-A x )K(x,y) = F( ia > x>+ 2 ^y > ). 
Particularly, if F G Mat s (K) and K G Mat s (A r ), then 
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(9) (I -A x )K 2 ( X ,y)= a / °° F( (^+C>+(ia,y> )p( (a,s>+(ia, ! r+C> )^ 

Acting on both sides of Equation (1) with the operator L x leads to the 
relation: 

(10) (a x - 1 a v )K(x,y) = * 2 a x a J °° F( ^±n>±i^> )K 2 (x, x + v )d V . 
Now we deduce an expression for a x K 2 (x, x + rf). From the definition of K 2 
the identity 

(11) (a x + 1 a z )K 2 (x,z) = (V, + \a z ) a f °° F( {a ' x+c> 2 + ^ z> )K(x,x + ()d( 

follows. Using (3.1) we get that 

(12) (a x + 1 a 2 )F(^±C^£>) = {(Tx+ A )F(»>) = ^F(QL (>..+<+.> , 
since fyo,- = fej and bj ±aj = bj for each j = 0, 2 r — 1. Therefore, in accor- 
dance with Proposition 4 the equality 

(13) O^ + ^K^) = (V x - V 2 ) a/ -F(^±^)K( a ;,a; + C)rfC- 
2F(^+£>)K(x,x) 

follows, where aiK(x,z) := a x K(x, z) and a 2 K(x,z) := a z K(x,z). We seek a 
solution K(x,y) satisfying the condition: 

(14) a y K(x,y) = i<r v K(x,y). 

From Equations (10, 13, 14) we get, that 

(15) (a x + a z )K 2 (x,z) = -2F(^¥>)K(x,x) 

+f V x ./ OO F(^f±^){ g J co [F( (6 ' 2x y K> )K 2 (x,x + t?)]^R, conse- 
quently, 

(16) (/-A x )[(<7 x +OK 2 (x,z)] = -2F(^>)K(x,x) = -2[(/-^)K(a;^)]K(x,x). 
If either r = 2 or F e Mat s (R) and K e Mat s (^4 r ), then one obtains 

from (3.3, 5, 6, 13) analogously the identities: 

(17) K- lC gK(x, „) = f (V- -V 2 ) (j F(^±^){ „ /» [F( ^± 2 ±C> ) K(x, *+ 
()}d(-2F(&^)K(x,x)}dr l 

= A x [(a x - a y )K(x,y)] - § F(^^)[F(^>)K(x, x)]^. 
Therefore, 

(18) (/ - A,)[K - <7 w )K(a;, y)] = -f [(/ - ^)K 2 (x, y)]K(z, x) + f , J °° < 
F( ^±p> ), F(^ft!Z>), K(x,x)>cfy. 

In the associative quaternion case ^4 2 = H or when F e Mat s (R) and K e 
Mat s (v4. r ) the last additive in Formula (18) vanishes, i.e. 

(18.1) (/ - A x )[(a x - a y )K(x,y)} = -\ [(/ - A x )K 2 (x, y)]K(x, x). 
On the other hand, one has the identities: 

(19) K 2 (x, x + Q = t ,JS° F( {b ' 2x T C> )K(x, x + £R and 
(19.1) (/ - A x )K 2 (x, z) = c j~ F(3z±p*>)K(x, x + CR 

-fa Jo" F(^^){F(^±^ )[ct/o oo F (^±^>)K(x,x+0^]dC}^ = 
<j Jo 00 F(^±f±^)[(7 - A x )K(x, x + CM 
= g / ° o F( (M+ 2 2+c> )F(^f^)dC, that is 
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(20) (I-A x )K 2 (x,z)= CT /~F(^±^)F(^C>K. 

If either r = 2, i.e. A 2 = H, or when F e Mat s (R) and K e Mat s (A), 
s G N for r = 2, s = 1 for r = 3, Equations (16, 18.1) imply that the 
functions K and K 2 are solutions of the system of partial differential equations: 

(21) (a x + a z )K 2 (x,z) = -2K(x, z)K(x, x) and 

(22) (a x - a z )K(x,z) = -%K 2 (x,z)K(x,x). 

The action on both sides of Equation (1) with the operator a x + \o y leads 
to the identities: 

(23) (a x + ia y )K(x,y) = a z F(z)\ z=(btX+y>/2 

+1 3 ^ ,SS° a Jo°° {F(^f±^)[F( ^ + 2 ^> )K(x, x + 0]d(}dv 
_E F( (te>) CT/ -[F(^fC>)K(a;, a ; + C)]rfC, 
since 

W y + 2 a x + K){F(^^)[F(^f ± ^)^(x,C)]} 
= 2a T; { F(^±^> ) [F( ( fc - 2 ^> C)] } 
for each A r valued function g(x,Q, consequently, 

(24) (I-A x )[(a x + i<r y )K(x,y)] = a z F(z)\ z={b , x+y>/2 -^[{I-A x )K(x, y)]K 2 (x, x). 
Then we use the condition 

(25) L 2 j F = for each j = 0, 2 r - 1, 
where 



(26) l, 2 .j = 2 ipjd tj + Vjihaj + 3 \i: ,<),,/:. + 3^ o'p),. + L o 



and act on both sides of Equation (1) with the operator 

L 2 = 2 <J t + al + 3 icr y a 2 + 3 \cr 2 a x + !(X 3 with = 0, i^o = 0, that gives 

(27) L 2 K(x, y) = E( 2 <r t + ^ + 3 ^ + 3 ^ + x^ 3 ) 
ct/ oo ^ f oo {F( ^±y> )[f{ (^> )K(a ;, g + C)]dC}^- 

To simplify appearing formulas we use the identity 

(28) (%+ 2 L 2 ) { F(^±fy>) [ F ( feg^±zz> ) g(a;> 0] } = 

V^f 1 ^ 2 <7 X + 2 <7 X Vj+3ja 2 V^^ 1 ^ 2 a x + 2 a x 1 <j x + V 2 )} 
{F( (^ )[F( (^> )g( ^ )]} 

= 3o- r ,( 2 <7 a 1 o- <: + 1 a c 2 cr z ){F(()\ (b, x +r, +y > \F(z)\ (b,2 X +c+v> g(x,Q}} 

^2 2 

for every A r valued function g(x,£), since the function F satisfies Conditions 
(25) and 

(28.1) a x l( x 2 y [f(x, y)g(x, y)\ = £a x + 2 a x ) [(\a y 2 a y + \a y \<j y )+(\<j 2 + 2 a 2 )] [f(x, y)g(x, y)\ 

and 

(28.2) a x [f(x, y)g(x, y)] = ?a x + 2 a x ) [( V x 2 a x+ 2 a x V X )+(V 2 + V 2 )] [f(x, y)g(x, y)] 
for each A r valued twice differentiable functions / and g. 

Then we have also from (3, 5): 

(29) (a x -a y )K 2 (x,y) = ( l a x - V,+ 2 a x ) a tf[F(Sb*±pt>)K(x, x+CM 
= 2 a xa f™[F(&^)K(x,x + ()]<!(, 

since d Xj F(x,y) = d yj F(x,y) for each x, y and j. 
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Hence, the identity is satisfied: 

(30) a x 2 a x [ J»F(^>)K(,,i + CK] = "2 2 a x [K(x,x)K(x,x)} - 
2( 1 a x - 1 a y )[K(x,y)K(x, x)}\ y=x - 2[a x , a y }K 2 (x,y)\ y=x , 

since [(a x + a y ), (a x - a y )\ = -2[a x , a y ] and 

o x V,[ CT / ~F(^f^)K(x,x + C)rfC] = ^{[K-^)K 2 (^y)]| 2/ =,} 

= (<7 X + (Tj,)^ - 0-^)^(2;,?/) I y=x 

= (cr x - <T y )(cr x + <r y )K 2 (x,y)\ y=x - 2[a x ,a y ]K 2 (x,y)\ y=x 

= (cr x - a y )[-2K(x,y)K(x,x)]\ y=x - 2[a x , a y ]K 2 (x, y)\ y=x 
due to Identity (16). On the other hand, 

[a x ,a y ]K 2 (x,y) = (V, V,- \ l a x + 2 a x V,- \ 2 a x ) a /-[F(^K±I>) K (x, x+ 
OR 

= ( 2 o x V,- V, 2 a x ) a J^[F(^±^)K(x,x + 0]dC, 
since ( 1 a :c — 1 a y 1 a x )[F( ^ b ' x+ ^ +y> )K(x, x + ()] = 0, consequently, 

[a^K^y)^ = i( V. V,- V s 2 a x ) a f™[F(&*^>)K(x,x + ()}d( 

and 

(31) (<7 a 2 a x + 2 a x a x ) [ a / °° F(^>)K(x, x+()d(] = -4 V, [K(x, x)K(x, x)} - 
4:( 1 a x - 1 a y )[K(x,y)K(x,x)]\ y=x . 

Using Equations (28 — 28.2) we rewrite Equation (27) in accordance with 
4(9 - 11) in the form: 

(32) (I-A x )[L 2 K(x,y)] = 

f(L 2 - 3 L 2 )„f °° (j Jo °{F(^f^)[F( ^ 2x+ 2 ^> )K(x, x + 0]dv}d( 

+f [{ l a x 3 a x + 3 a x l a x + 2 a x 3 a x + 3 a x 2 a x )\a y + \a 2 3 o x +( 1 a x 2 a x + 2 a x 1 a x + 1 a 2 x + 2 a 2 x ) 3 a x ] 

u/ oo J ^{F(g»)[F( ^ + ''> )K(x,. + C)]^K 

+f[(V.+ V* + S) 3 ^( CT / °° ct/o oo {f( (m^ )[f( (m^c + ,> ) 

K(x,x + C)]d V }dC 

= ~ f t 1 ^ ^ + ^ V] - fo°° F(^) l 2= (Md^> [F(0 1 5= (m^c> K(x, x + CM 

-f 3 ^./o ,w F(fi^>)[F(a^)K(x,x + C)]dC 
+3p 3 ^ CT J- g ; «»{F(^m> ) [F( )K(x, x + C)]^R 
since 

(^ + is 2 )[ F ( (M T" > ) F ( (fc ' 2a:+ 2 C+T;> )] = ^( ib ' x+ r y> )H ib ' 2x+ 2 +v> )]- 

If F e Mat s (R) and K e Mat s (A r ), 2 < r < 3, s E N for r = 2, s = 1 
for r = 3, then using Equation (11) we deduce that 

(33) (I- A x )[L 2 K(x,y)} = -f (V, 2 a x + V, V,)^)^^, *)] - 
f 3 ^(V K + 3 a x ){F(&±v>) ff /nF(^)K(^ + ()K}. 2 

From Equations (16,31,33) the identity 

(34) (/_^ a: )[L 2 K(x, 2 /)] = {^[-f(/-A x )( ( 7 x + ia y )K(x,y) 
- 3 4[(I - A x )K(x,y)]K 2 (x,x)]K 2 (Z,Z)}\s =x 

follows. Expressing F from Equality (1): 
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F(^y>) = F(x,y) = (I-A x )K(x,y) 
and using Formulas (21) and 4(7,8) and the invertibility of the operator 
(I — A x ) and 7(13) (see Remark 3.2 also) we get 

(35) L 2 K(x,y) = {[-f (a x + ia y )K(x, y) — %£-K(x, y)K 2 (x, x)](-K 2 (x, x)} 
-f{K(xj)[(-2V)(K(x,x)K(x ) x))]-2K(x,y)[(V,-| S )(K( a : l!/ )K(x ia; ))]|, =x } 
= 3p(V a + l( T y )K(x,y)K 2 (x,x) + ^[K(x, y)K 2 (x, x)]K 2 (x, x) 

+3pK(x, y) [ 2 a x (K(x, x)K(x, x))]+3pK(x, y) [{ x a x - \<J y )(K(x, y)K(x, x))]\ y=x . 
Using Formula (22) permits to simplify Formula (35) in the considered 
case to: 

(36) L 2 K(x,y) = 3p[(a x +a y )K(x,y)}K 2 (x,x)+3pK(x,y)[ 2 a x (K(x,x)K(x,x))}. 
Putting g(x, t) = K(x, x, t) one gets particularly for a — \<j from Equation 

(36): 

(37) ( 2 a t + a 3 x )g(x, t) = 3p[a x g{x, t)]g 2 {x, t) + 3pg{x, t) [ 2 a x (g(x, t)g{x, t))}. 
The latter two partial differential equations are non-linear even over the 

real field. Particularly, if 2 a t = d/dt and a x = d/dx Equation (37) gives 
the modified Korteweg-de-Vries equation: 

dg(x,t)/dt + d 3 g(x,t)/dx 3 — 6pg 2 (x,t)dg(x,t)/dx = with t — t and 
x = x . 

Find particular solutions of (36), when 

2 a t f = (df/dt ) and a x g(x) = E^i i*[dg(x)/dx j ], using 

F(^>),0 := (nQ-WfiU+l+Tho+l, -^n+no+Vn+no) with £ = x{nt)~ l l 3 

and n = y(nt) -1 / 3 , — 1 < n , 2 < n, n + n < 3, t = t . While with the 
help of the operators 

L 2J g(x, y, t) = li*d t g(x, y, t) + (d Xj a 2 x + d Vj a 2 + 3d Xj a 2 + 3d yj a 2 x )g{x, y, t) 
and the conditions L 2 jF = for each j — no + 1, no + n one obtains the 
partial differential equations 

/ + &/ + 0&£ffl%/ = o 

for each j = n + 1, ....,n + n , where / = /(£„ +i, ...,£ no+n ), d Xj = d/d Xj . 
Solving these partial differential equations we get the solution K(x,y) = 
(nO~ 1/3 K(Cn +i, -,fno+n,»to, -,»7no+n) from the integral equation 

where X (£, 77) = «(£„„+!, £n +n, ?7o, Vn +n), £ = + - + &»3- 

Each solution of the modified Korteweg-de-Vries equation generates a so- 
lution of Kadomtzev-Petviashvili partial differential equation due to Miura's 
transform: g = —v 2 — a x v. In the considered non-commutative case it takes 
the form: 

2 <T t (v 2 ) + 2 a t a x v + (a x (v 2 ) + o%6) - 3p(a x (v 2 ) + a 2 x v) (v 4 + v 2 a x v + (a x v) v 2 + 
(vxv) 2 ) - 3p(v 2 + a x v)[ 2 a x ((v 2 + a x v)(v 2 + a x v))}. 

Thus we have demonstrated the following statement. 
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9.1 Theorem. Partial differential Equation (36) with ip = iip = over 
the Cayley-Dickson algebra A r , 2 < r < 3, has a solution given by Formulas 
(1 — 3, 25, 26), when the appearing integrals uniformly converge by parameters 
as in Proposition 4 and the operator (I — A x ) is invertible and F G Mat s (R) 
and K G Mat s (A r ) with s G N for r = 2 and s = 1 /or r = 3. 

10. Example. Now we consider the pair of partial differential operators 

(1) L\ = a x + o~ y and 

(2) L 2 = i<r t + <?l + ua y a x + a 2 , 

where u G R is a real constant, and the integral equation: 

(3) K(x, y) = F(x, y) + p , £° F(z, y)K(z, z)dz. 

Acting on the both sides of the latter equation with the operator L 2 and 
using the condition 

(5) L 2 jF(x : y) = for each j, where 

(6) L 2J F(x,y) = [iTpjd tj + 

uJ2s, P : i s i p=ij ^s^p{d ys d Xp + {i p i s )(i s i p )d yp d Xs ) + tt,- o (^ + <j 2 )]F(x,y), 
where 7T,- : X — >■ X, is the R linear projection operator (see §1.2). We infer 
that 

(7) L 2 K(x,y) = pHo-t + al + uo-y^ + o-l + jo-t) a f™ F{z, y)K(x, z)dz 
= p( 2 a t + o-l + u 1 o y o x - 1 o- 2 z -u 1 a y 1 a z ) a J™ F(z, y)K(x, z)dz, 

since (i<7 t + a 2 )F(z, y) = — {o- 2 + ua y a z )F(z, y). In view of Proposition 4 and 
Corollary 5 we get the formula: 

(8) L 2 K(x,y)=h + I 2 , 
where 

(9) h = p( \°t + o~ 2 x + u l a y a x ) a J™ F(z, y)K(x, z)dz 
= p( 2 a 2 x + u x Gy 2 a x ) (j F(z, y)K(x, z)dz 

-pa x [F(x,y)K(x,x)]-p 2 a x [F(x,y)K(x,z)}\ z=x -upa y [F(x,y)K(x,x)] and 

(10) -I 2 = pCal + u'ay V 2 ) „£° F(z,y)K(x,z)(fe 
= P( 2 ^~ M 2 °z) aI^°F(z,y)K(x,z)dz 

-p 1 a x [F(x,y)K(x,x)] + p 2 a z [F(x,y)K(x, z)}\ z=x - upa y [F(x,y)K(x,x)]. 
Then using the second condition 

(11) L 1 F(x,y) = 0, i.e. F(x, y) = F(^p>), 

it is possible to simplify certain terms appeared above: 

(12) <j y 2 <j xa f™F(z,y)K(x,z)dz 

= 2 o- z 2 a x a J™F(z,y)K(x,z)dz+ 2 (T x [F(x, y)K(x, z)}\ z=x and 

(13) \ 2 a z a f™F(z,y)K(x,z)dz = - V, 2 a z c F(z, y)K(x, z)dz 
= 2 a 2 „f™F(z,y)K(x,z)dz + 2 a z [F(x, y)K(x, z)\ \ z=x . 
Therefore, Equations (7 — 13) lead to the identity: 

(U)L 2 K(x,y) = p( 2 a t + 2 a 2 x +u 2 a z 2 a x +(u-l) 2 a 2 z ) a F(z, y)K(x, z)dz 
-(u+2)p 2 a x [F(x,y)K(x,x)]+up 2 a x [F(x,y)K(x, z)]\ z=x +up 2 a z [F(x,y)K(x, z)]\ z=x . 
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Take u = 2. If F(x, y) G Mat s (H) and K(x,y) G Mat s (A r ) for each 

x, y G £7, 2 < r < 3, with s G N for r = 2 and s = 1 for r = 3, then Formulas 
(3, 14) and 2(7) and 7(2, 13) and 5(7) for the invertible operator (/ — A x ) 
imply, that 

(15) k<7 t + ul + 2 2 a y 2 a x + a 2 y )K(x,y) = —2pK(x, y) [a x K(x, x)] . 
Putting g(x,t) = K(x,x) we get particularly on the diagonal x = y: 

(16) (i(j t + al)g(x,t) = -2pg(x,t)[a x g(x,t)}. 

If u G X r := {z G A r : Re{z) = 0} and = 0, then 2 <T x (u(a;, t)u(x, t)) + 
u(x,t)(a x u(x,t)) = J2j>o u j( x ,t)(du(x,t)/dxj)ipj- 
This approach has the following application. 

A non- isothermal flow of a non-compressible Newtonian liquid with a 
dissipative heating is described by the system of partial differential equations: 

(17) pE^ujdu/dx, = -\n°lu + 

(18) div(u) = Ej 3 =i duj/dxj = 

(19) pc p E| =1 UjdT/dXj = -\a 2 x T + 2/iJ 2 , 

when a density p and a dynamical viscosity /z are independent of coordinates 

xi, X2, X3, where 

h = \ Elk^iduk/dxj + duj/dxkf 
denotes the tensor of deformation velocities, T is the temperature, also A 
and c p are physical constants. Here we take u = u\i\ + U2I2 + "^323 with real- 
valued functions mi,m 2 ,U3 and a x f(x) = Y?j=ii){.df(x)/dxj) ) x = xoio + 
X\i\ + X2H + ••• + Xyr-iiw-i G A r , r = 2 or r = 3, x and Xj with j > 3 can 
be taken constant, for example, x = 0, Xj G R for each j, \<j t = iV'o an( i 
choose suitable constants iipo and p. 

For usual liquids like water the dependence of the dynamical viscosity 
on the temperature is described as fi = /io(T /T) m , where m > 0. For very 
viscous Newtonian liquids like glycerin the function /i = /i exp[— (3(T — T )] is 
usually taken with empirical constants /io, T and /3 (see also |18j). Using the 
method described above it is possible to find solutions satisfying Equation 
(20) and take into account Condition (21) with the help of F(x, x) and then 
K(x, x) satisfying it. 

Put g = v + u, where u G X 2 := {z G A2 '■ Re(z) = 0} and v G A r X 2 , 
where 2 < r < 3. If g is a solution of partial differential equation (16) and 

v(au) + u(cxv ) = 2 a(uu + vv), 
then m is a solution of partial differential equation (17). Condition (18) means 
that Re(au) = or au + (cru)* = 0. Calculating J 2 and solving Equation 
(19) one also calculates the temperature. 

The result of this section can be formulates as the following theorem. 

10.1 Theorem. Partial differential Equation (15) over the Cayley-Dickson 
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algebra A r with 2 < r < 3 has a solution given by Formulas (1,3,5,6, 11), 
when the appearing integrals uniformly converge by parameters on each com- 
pact sub-domain (see Proposition 4) and the operator (I — A x ) is invertible 
and F E Mat s (R) and K E Mat s (A r ), s G N for r = 2, s = 1 for r = 3. 

11. Remark. The results presented above show that the method of 
commutators of integral operators with partial differential operators over 
complex numbers becomes more powerful with the use of Cayley-Dickson 
algebras and the non-commutative line integral over them and permits to 
solve more general non-linear partial differential equations. It is planned to 
develop further this method for partial differential equations with variable 
coefficients and also for generalized functions and equations with generalized 
coefficients. 
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